High Harmonic Generation via Continuum Wave-Packet Interference 
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High-order harmonic generation (HHG) is investigated theoreticahy in the over-the-barrier ion- 
ization (OBI) regime reveahng the strong signature of interference between two separately ionized 
and separately propagating free wave packets of a single electron. The interference leads to the 
emission of coherent light at a photon energy corresponding to the kinetic-energy difference of the 
two recolliding electron quantum paths, thus complementary to the well-known classical three- step 
picture of HHG. As will be shown by time-frequency analysis of the emitted radiation, the process 
entirely dominates the coherent HHG emission after the atomic ground state has been depleted by a 
strong field. Moreover, it can be isolated from the continuum-bound harmonics via phase-matching. 

PACS numbers: 42.65.Ky, 42.79.Nv 



High-order harmonic generation (HHG) is a key pro- 
cess in ultrafast science. As a fundamental example of 
strong- field electron dynamics, it is fascinating both from 
a fundamental physics perspective — coherent recolli- 
sion of sub-femtosecond ionized electronic wavefunction 
within a fraction of an optical cycle of the intense driv- 
ing laser [THl] — but also from a technological point of 
view enabling the creation of attosecond-duration flashes 
of light [5, 6 for the time-resolved study of the fastest 
processes in physics and chemistry (for current reviews 
see [II IBj). HHG is currently well- understood in terms 
of the three-step model [U [3] : An electronic wave packet 
is partially ionized, accelerated in the laser field, and re- 
turns to its parent ion. At that point, the ionized and 
bound-state portions of the electronic wave packet in- 
terfere giving rise to a strong, coherent high-frequency 
dipole response that can lead to the emission of an HHG 
photon along with the recombination of the electron into 
the bound state. The coherent HHG emission from this 
continuum-bound (CB) mechanism is practically absent 
as soon as the atom is fully ionized by the laser pulse [4] . 
Then, recombination to the ground state can only pro- 
ceed by incoherent spontaneous emission of radiation 
that scales linearly with the number of atoms N in the 
gas target rather than ex N'^ as in the case of the very 
efficient typically considered coherent process. 

However, the impact of free portions of the electronic 
wave packet alone has not been fully recognized yet. 
Continuum-continuum (CC) transitions in atomic HHG 
have been studied with focus on the Bremsstrahlung 
emission of single wave packets recolliding with a bare 
core |9] and the interaction with a continuum wave packet 
ionized just shortly before the recollision time [lOj lllj . 
Also, in a recent theoretical work, it was shown that rec- 
ollisions of two free electrons after strong-field double- 
ionization can lead to weak HHG emission beyond the 
usual cutoff 12 . 



of atomic HHG by predicting the efficient production of 
coherent soft-x-ray light from the interference of spec- 
trally distinct free single- electron wave packets. It is 
shown both analytically and numerically that even for 
the case of a fully depleted single-electron atomic sys- 
tem (no bound states populated) a characteristic coher- 
ent polarization response can be obtained correspond- 
ing to the energy difference of two distinct continuum 
electronic states. Investigating the temporal structure of 
emission reveals that the process dominates the tradition- 
ally known CB interference of the electron recombining 
with the bound states at times after the ground-state 
population has been depleted. 

In most general terms, the continuum-continuum (CC) 
interference mechanism presented here provides a fun- 
damental physical picture: wave function splitting or 
spreading, subsequent simultaneous recollision with dif- 
ferent energies, and core-mediated transition with pho- 
toemission at the difference energy. This approach thus 
generalizes the traditional picture of CB HHG to include 
CC transitions. 

We start out by describing the key idea considering the 
case of an electron-atom collision. The coherent part of 
the emitted radiation can be calculated via the expecta- 
tion value of the acceleration using Ehrenfest's theorem: 



a{t) = -{*(i)|VF|*(t)) 



(1) 



with the ionic potential V and atomic units (a.u.) being 
used throughout. Note that this approach, in consistence 
with earlier and current atomic HHG theory, considers 
only the coherent part of the emission response and ne- 
glects incoherent spontaneous emission which cannot be 
phase-matched. We briefly describe the general condi- 
tions required for the emission of radiation: If |^(t)) is 
an eigenstate of the Hamiltonian, the expectation value 
is time-independent. Harmonic emission can only oc- 
cur for a linear combination of at least two eigenstates: 
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dependent part a{t) = —aia2{ei\\/V\e2)e' 
shows the resulting dipole oscillates with the difference 
energy between the two states. CB HHG can be un- 
derstood in this context: A recolliding electronic plane 
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FIG. 1: (color online) Electron-atom collision and analy- 
sis of the dipole acceleration response. Windowed Fourier 
transforms on a logarithmic scale are shown for four cases: 
a) monochromatic electron wave packet (47.8 a.u. kinetic en- 
ergy) colliding with a partially ionized atom, b) same as a) 
but for a fully ionized atom, c) a bichromatic electron wave 
packet (47.8 a.u., 33.1 a.u. kinetic energies) colliding with a 
fully ionized atom, d) same as c) but for a partially ionized 
atom. In c) and d) the signature of cc electron interference is 
observed at the difference kinetic energy of 14.7 a.u. 



wave (energy ei) interferes with the bound wave packet 
(energy —Ip) and the difference energy (ei + /p) is emit- 
ted. Moreover, if a bichromatic electron wave sweeps 
over the potential, it would lead to light emission at the 
kinetic energy difference. This intuitive model strictly 
holds only for time-independent Hamiltonians but its pic- 
ture remains in atomic HHG with collision times typically 
short on the time scale of a laser period. 

To numerically illustrate this difference-frequency 
mechanism, we solve the time-dependent Schrodinger 
equation for a 1-dimensional hydrogen a tom us ing a 
smoothed Coulomb potential V{x) = 1/Vx^ + a^ with 
the smoothing parameter a = 1.4039 chosen to match 
the hydrogen ionization potential. We consider 4 cases: 
a) a monochromatic free-electron wave packet colliding 
with a partially ionized atom (the conventional HHG sce- 
nario), b) the same as a) but for a fully ionized atom, 
c) a bichromatic free electron colliding with a fully ion- 
ized atom, and d) a bichromatic free-electron wave packet 
colliding with a partially ionized atom. We analyze the 
dipole acceleration expectation value by means of a win- 
dowed Fourier transform (Gabor transform) with a tem- 
poral window size of 3 a.u. The results are shown in 
Fig.[l] 

In case a), we observe the dipole response at a fre- 
quency corresponding to the sum of electron kinetic en- 
ergy and ionization potential of the atom £^kin + ^p, well- 
known from traditional HHG. In case b), no high-energy 
dipole response can be observed. The low energy signal 
corresponds to the Bremsstrahlung spectrum [9 arising 



due to its non- vanishing spectral width and in fact can be 
understood as a special case of the mechanism discussed 
here. Case c) exhibits a dipole response at a slightly 
lower frequency than the one obtained in a), while case 
d) shows a total of 3 different frequencies. The absence of 
any high-energy emission in case b) arises from the lack 
of interference as only one electronic state takes place in 
the interaction in consistence with an earlier theoretical 
study [4 . Emission in c) is created by the interference 
of two free-electron wave packets with kinetic energies 
of £^kin,i and ^kin,2 that create a time-dependent dipole 
moment corresponding to a beating frequency of the dif- 
ference ^kin,i — ^kin,2 inside the non-linear potential set 
up by the Coulomb field of the atom. This is the origin 
of the coherent CC electron emission process for HHG 
described in more detail below. Case d) shows emission 
due to both free-electron wave packets, each interfering 
with the bound state and amongst themselves, giving rise 
to three distinct dipole oscillation frequencies. Note that 
CC transitions also occur when the two wave packets 
meet each other in the potential with opposite (differ- 
ent absolute) momenta but with a lower probability as a 
consequence of the higher momentum transfer required. 

In the following, it is shown that HHG radiation con- 
tains a signature of CC electron wave packet interference 
when approaching the over-the-barrier (OBI) regime. As 
a more realistic model, in this case, we use a hydrogen 
atom with V(r) = — 1/r and solve the 3-dimensional 
time-dependent Schrodinger equation in a strong field 
but in principle, any other atom, ion or molecule could 
be chosen. The peak intensity of lO^^W/cm^ of the 
laser pulse (shown in Fig. |2^) is chosen such that al- 
most complete depletion of the ground state occurs on 
the leading edge of the pulse. The evolution of the 
ground state population is shown in Fig. [2|i) via no{t) = 
exp (— /_ drw{r)) where w{r) is the ionization rate be- 
ing a combination of the empirical formulas in Ref. [13] . 
This way, CC transitions become the dominating HHG 
mechanism for large parts of the pulse. The laser in- 
tensity and recollision energies are higher than in most 
of the current HHG experiments but our results are not 
specific to these intensities. For solving the TDSE, the 
code QPROP [14 is used which expands the wave func- 
tion in spherical harmonics and propagates them within 
a Crank-Nicolson scheme. As simulation parameters, 
we used a time step of At = 0.01 a.u., a radial grid of 
^max = 800 a.u. with a spacing of Ar = 0.1 a.u. and 
a highest angular momentum of L = 140. Again, we 
calculate the windowed Fourier transform of the dipole- 
acceleration response shown in Fig. [2b). 

To analyze the individual time-resolved frequency com- 
ponents contained in the dipole response, we use New- 
ton's equations to calculate classical electron trajectories 
in the laser field originating from the ion. In Fig. [2]3, 
dashed black line, we show the classical continuum- 
bound (recombination) energy of £^kin,i,2 + ^p as well 
as the energy difference £^kin,i — ^kin,2 of the two re- 
turning classical trajectories (solid red), each vs. the 
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FIG. 2: (color online) Time- frequency analysis of HHG show- 
ing the signature of CC wave-packet interference, a) Driver 
pulse used in the simulation (solid line) and the reliable part 
of the ground-state population (dashed line) calculated via 
[13| . b) Windowed Fourier transform of the HHG emission 
on a logarithmic scale. Dashed black lines: classically calcu- 
lated kinetic energies of electrons returning to the ion. Solid 
red line: difference between the two black dashed curves. 



ollision process and harmonic emission by an analytical 
theory framework based on the strong-field approxima- 
tion (SFA). The standard approach to HHG within the 
SFA where the dipole moment is calculated [2 is not suf- 
ficient since two plane continuum waves with different 
momenta have a vanishing dipole moment. As in the 
SFA no distortion of the plane waves by the potential of 
the atom is included, the impact of the potential is im- 
posed by directly calculating the nonlinear dipole accel- 
eration responsible for HHG emission via the expectation 
value of the acceleration (Eq. IT]). This way, the interac- 
tion of the electron with the core is included which is 
needed for the momentum conservation of the CC tran- 
sition [15 . Within the dipole approximation (DA), the 
electronic wavefunction can be written as 



l*(t)) = \Mt)) + \4>c{t)) 

= Vm(t)\0)e"-'-i fd^p\p + A{t)) 



(2) 
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time of recollision. We compare those to the dipole- 
response obtained from our quantum calculation: Until 
time t = 150 a.u., HHG emission is caused by CB transi- 
tions that are still rather strong because the bound state 
has not been fully depleted (compare Fig. [2] a). After 
time t = 150 a.u., we can clearly observe both the signa- 
ture of ground-state recombination (CB) as well as CC 
interference by comparing to the classically calculated 
return energies. Interestingly, the CC component of the 
dipole response is the dominant contribution for several 
half-cycles between t = 150 a.u. and t = 320 a.u. This 
can be understood from the fact that depletion of the 
ground state occurs around t = 150 a.u. and thus, inter- 
ference with and recombination to the ground state can 
no longer take place. Instead, the two parts of the wave 
packet which have been ionized within the first two half 
cycles steadily interfere with each other in the atomic 
potential for the remainder of the laser pulse. Note that 
the electronic interference described here is different from 
the optical interference of the harmonic field caused by 
the well-known short and long trajectories. 

The significance of the signal is pointed out by estimat- 
ing the photon yield in a small energy window between 
12 a.u. and 18 a.u. Assuming a phase-matched gas vol- 
ume with radius 100 /im and length 1 mm at a density 
of 10^^ /cm^ and emitting into a solid angle of 10~^, we 
calculate a yield of 10^ photons per shot within the time 
window of t = 110 a.u. and t = 155 a.u. originating from 
CB transitions. Between t = 155 a.u. and t = 200 a.u., a 
photon number of 10~^ is emitted within this bandwidth 
mainly caused by CC transitions. The typical kHz rep- 
etition rates, therefore, allow for a measurement of this 
effect that could potentially be enhanced by using larger 
atoms/molecules with stronger core potentials. 

In the last part of this letter, we describe the CC rec- 



where S{p,t,t') = jLdT{p + A(r))^/2 is the classi- 
cal action. Inserting pi) into (IT]) yields two contribut- 
ing terms [11 : The CB transition can be identified as 
acsit) = -{(l)B{t)\S/V\^c{t)) -^c.c. and the CC transi- 
tion as acc{t) = -((^cWlVVl^cW). 

The saddle-point approximation (valid for uj <^ Ip ^ 
Up) is successively applied to both expressions. This 
way, the momentum integration is carried out first with 
saddle points defined via \/pS = yielding ps{t^t^) = 
— J^, drA{r)/{t — t') for a field polarized linearly along 
the X-axis. The time integration is more involved as the 
ionization matrix elements exhibit a singularity at the 
saddle points of the t' (and t") integration defined by 
dt'S{ps{t^t^)^t^t^) — Ip = but can be treated in a simi- 
lar way as in Ref. [16 . Finally, we obtain 



acB{t) 



acc{t) 



= -Re ^^^P^ Y^ ^no(t)no(t')wcorr(t')D(t, t') 

t' 
X (0|5^F|p,(t,tO + A(t))e-^(^(^^(^'^')'^'^')+^-(^-^')) 
(8/p)^/2 
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-i{S{psit,t"),t,t")-S{ps{t,t'),t,t')+Ip{t'-t")) 



where D{t,f) = ^{2mf/{t - t'f /{d^ /dt'^S{ps{t, t'),t, t')), 
Wcorr{t) = \/u){t)/wK{t) and t' < t" . The ionization 
correction factor Wcorr{t) was inserted ^ as the 
non-Coulomb corrected SFA contains only the Keldysh 
ionization rate wxit) which is not sufficient in this 
regime. The matrix elements {p\dxV\p') and {0\dxV\p) 
can be calculated analytically for a spherically sym- 
metric Coulomb potential V{r) = —y/ajr with the 
ground state wave function V^is = a^/^/y^e~^^ where 
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FIG. 3: (color online) CC HHG by using an analytical ap- 
proach based on the SFA for the same field as in Fig. [2] (a) 
is the windowed Fourier transform of the acceleration within 
a chosen time window in excellent agreement with the TDSE 
results of Fig[2] (b) shows single-atom spectral photon yield 
per solid angle per shot of the same time window as in (a). 



_ i^/a 



and 



a = 2Ip which yields {Px\ - da^Vlq^,) - ^^2 p^_q^ 

{Px\ - d,V\0) = iV2a^/^{p, - arctan(p,/a))/(^pf). 

The analytically calculated dipole acceleration is 
shown in Fig. |3J for the same electric field and in the 
same manner as for Fig. [2] The excellent agreement 
with the TDSE results shows that the process is cor- 
rectly described by the SFA model. The SFA model 
allows to extract the phase of the emitted harmonics 
cj) = S{p,{t,t"),t,t")-S{p,{t,t'),t,t')+Ip{t' -t")-ujHt 
converted to (/)(/) = acci as a function of pulse peak in- 
tensity /, which reveals a striking difference for CB as 
compared to CC transitions. For the time window con- 
sidered in Fig. [3J we find 6 contributions to the harmonic 
energy at ujh = 8 a.u.: Two CB transitions from the 
wave packet ionized at the first peak of the laser pulse 
(<^CB,i,s = +0.6 X 10"^^cm^/W and acB,i,i = +1-5 x 
10~^^cm^/W for the short and long trajectory, respec- 
tively), two CB transitions from ionization at the second 
peak {acB,2,s -^ and acB,2,i = +1-4 x 10"^^cm^/W) 



and the two CC transitions {acc,s = — 0.6x 10~-^^cm^/W 
and acc,i = -0.1 x 10~^^cm^/W). As a result, and as 
can be understood by the difference-energy mechanism, 
the sign of a for CC vs. CB harmonics is found to be dif- 
ferent and thus allows to experimentally isolate the CC 
harmonics. 

The impact of our results could be manifold: First of 
all, CC HHG may qualitatively advance and complete 
tomographic molecular imaging: Instead of sensing the 
orbital shape of the active electron [18] , the effective po- 
tential could be measured. This is because the detailed 
structure of the potential enters the CC transition mainly 
via the matrix element {p'\\/V\p") which is the Fourier 
transform of the gradient of the effective potential, which 
is not restricted to an ionic potential. This could be par- 
ticularly interesting for larger molecules, where the entire 
molecular potential, i.e. mainly the positions of all nu- 
clei, would be mapped out, not just the few highest-lying 
molecular orbit als. 

Moreover, it is known that an electron ionized from one 
isolated atom cannot coherently recombine with another 
atom to produce CB HHG, as there is no phase relation 
among isolated atoms. However, the phase of the CC 
transition does not depend on any bound-state portion 
of a wavefunction, only on the shape of the potential. It 
would allow for coherent interferences in another atomic 
or more complex potential even if they are some fixed 
distance away. 

Also, CC transitions contribute to harmonic energies 
below Ip on which currently a frequent discussion exists 
about their underlying mechanism \T9\. 
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